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Abstract 



O 
O 

o 

r~| ' We study the semileptonic transitions — > r/c) J/i^^ D, D*, B, B*, Bg, B* in 

the framework of a relativistic constituent quark model. We use experimental 
' data on leptonic J/^p decay, lattice and QCD sum rule results on leptonic Be 



decay, and experimental data on radiative rjc transitions to adjust the quark 



' model parameters. We compute all form factors of the above semileptonic 

Sc-transitions and give predictions for various semileptonic B,, decay modes 
including their r-modes when they are kinematically accessible. The impli- 
cations of heavy quark symmetry for the semileptonic decays are discussed 
and are shown to be manifest in our explicit relativistic quark model calcu- 
lation. A comparison of our results with the results of other calculations is 
performed. 
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I. INTRODUCTION 



Recently, the observation of the bottom-charm Be meson at Fermilab Tevatron has been 
reported by the CDF Collaboration The Be mesons were found in the analysis of their 
semileptonic decays, Bf —>■ J/ipl^X. Values for the mass and the lifetime of the Be meson 
were given as M(5J = 6.40 ± 0.39 ± 0.13 GeV and t^B^) = 0.46l|]:}^(stat) ± 0.03(syst) ps, 
respectively. The branching fraction for B^ J /ip I v relative to that for Be —* J / ipK was 
found to be 

a{Be) X Bi{Be ^ J/i^lv) _ .32+0.041. ,,^ , . 03irsvstl+°-°32 
a{B) X Bi{Be ^ J/i^K) ~ U-l^^-o.osr^stat) ± 0.031(syst)_o.o2o • 

The study of the Be meson is of great interest due to some of its outstanding features. It 
is the lowest bound state of two heavy quarks (charm and bottom) with open (explicit) flavor 
that can be compared with the charmonium (cc-bound state) and the bottomium (66-bound 
state) which have hidden (implicit) flavor. The states with hidden flavor decay strongly and 
electromagnetically whereas the i?c-meson decays weakly since it is below the SD-threshold. 
Naively it might appear that the weak decays of the i?c-meson are similar to those of the 
B and D mesons. However, the situation is quite different. The new spin-flavor symmetry 
arises for the systems containing one heavy quark when the mass of the heavy quark goes 
to infinity [^. It gives some relations between the form factors of the physical processes. 
The deviations from heavy quark symmetry are large for the D meson and negligibly small 
for the B meson. On the contrary, in the case of the Be meson a consistent heavy quark 
effective theory (for both constituent quarks) cannot include the heavy flavor symmetry . 
However, the residual heavy quark spin symmetry can be used to reduce the number of 
independent semileptonic form factors at least near the zero recoil point [Q. 

In the naive spectator model, one would expect that T{B^ ^ r(i?) + T{D) which gives 
T{Be) ~ 0.3 ps, i.e. 1.5 times less than the central CDF value. The dominance of the c — >■ s 
transition will have to be investigated in future analysis when more data becomes available. 
Thus a reliable evaluation of the long distance contributions is very important for studying 
the weak Be decay properties. 

The theoretical status of the i?c-meson was reviewed in 0. In this paper we focus 
on its exclusive leptonic and semileptonic decays which are sensitive to the description 
of long distance effects and are free of further assumptions, as for example, factorization 
of amplitudes in non-leptonic processes. Our results on the semileptonic transition form 
factors can of course be used for a calculation of the nonleptonic decays of the i?c-meson 
using the factorization approach. 

The exclusive semileptonic and nonleptonic (assuming factorization) decays of the Be- 
meson were calculated before in a potential model approach Q]. The binding energy and 
the wave function of the Be-v[veson were computed by using a flavor-independent potential 
with the parameters fixed by the cc and hh spectra and decays. The same processes were 
also studied in the framework of the Bethe-Salpeter equation in [^, and, in the relativistic 
constituent quark model formulated on the light-front in [^. Three-point sum rules of QCD 
and NRQCD were analyzed in @,|l^ to obtain the form factors of the semileptonic decays 
of 5+ ^ J/iP{r]e)l+iy and B+ BsiB*)l+iy. 

As shown by the authors of 0, the form factors parameterizing the Be semileptonic 
matrix elements can be related to a smaller set of form factors if the decoupling of the 



2 



spin of the heavy quarks in Be and in the mesons produced in the semileptonic decays is 
exploited. The reduced form factors can be evaluated as overlap integral of the meson wave- 
functions obtained, for example, using a relativistic potential model. This was performed in 



111 , where the Be semileptonic form factors were computed and predictions for semileptonic 
and non-leptonic decay modes were given. 

In this paper we employ the Relativistic Constituent Quark Model (RCQM) [jl2[ for the 
description of Be semileptonic meson decays. The RCQM is based on an effective Lagrangian 
describing the coupling of hadrons H to their constituent quarks the coupling strength of 
which is determined by the compositeness condition = |T3] where is the wave 
function renormalization constant of the hadron H. Z]/^ is the matrix element between 
a physical particle state and the corresponding bare state. The compositeness condition 
Zh = enables us to represent a bound state by introducing a quasiparticle interacting 
with its constituents so that the renormalization factor is equal to zero. This does not mean 
that we can solve the QCD bound state equations but we are able to show that the condition 
Zh = provides an effective and self-consistent way to describe the coupling of the particle to 
its constituents. One starts with an effective Lagrangian written down in terms of quark and 
hadron variables. Then, by using Feynman rules, the S-matrix elements describing hadronic 
interactions are given in terms of a set of quark diagrams. In particular, the compositeness 
condition enables one to avoid a double counting of hadronic degrees of freedom. This 
approach is self-consistent and all calculations of physical observables are straightforward. 
There is a small set of model parameters: the values of the constituent quark masses and 
the scale parameters that define the size of the distribution of the constituent quarks inside 
a given hadron. This distribution can be related to the relevant Bethe-Salpeter amplitudes. 

The shapes of the vertex functions and the quark propagators can in principle be found 
from an analysis of the Bethe-Salpeter and Dyson- Schwinger equations, respectively, as 
done e.g. in |T^. The Dyson-Schwinger equation has been employed to entail a unified and 
uniformly accurate description of light- and heavy-meson observables [l^ . In this paper we, 
however, choose a more phenomenological approach were the vertex function is modelled 
by a Gaussian form, the size parameter of which is determined by a fit to the leptonic and 
radiative decays of the lowest lying charm and bottom mesons. For the quark propagators 
we use the local representation. 

The leptonic and semileptonic decays of the lower-lying pseudoscalar mesons (vr, K, 
D, Ds, B, Bs) have been described in Ref. Ijl^ in which a Gaussian form was used for 
the vertex function and free propagators were adopted for the constituent quarks. The 
adjustable parameters, the widths of Bethe-Salpeter amplitudes in momentum space and the 
constituent quark masses, were determined from a least square fit to available experimental 
data and some lattice determinations. We found that our results are in good agreement 
with experimental data and other approaches. It was also shown that the scaling relations 
resulting from the spin-flavor symmetries are reproduced by the model in the heavy quark 
limit. 

Using this approach we have elaborated the so-called Relativistic Three-Quark Model 
(RTQM) to study the properties of heavy baryons containing a single heavy quark (bottom 
or charm). For the heavy quarks we used propagators appropriate for the heavy quark 
limit. Physical observables for the semileptonic and nonleptonic decays as well as for the 
one-pion and one-photon transitions have been successfully described in this approach [jl8 
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Recently, the RTQM was extended to include the effects of finite quark masses [0. We 
mention that the authors of have developed a relativistic quark model approach to the 
description of meson transitions which has similarities to our approach. They also use an 
effective heavy meson Lagrangian to describe the couplings of mesons to quarks. They use, 
however, point-like meson-quark interactions. Loop momenta are explicitly cut off at around 
1 GeV in the approach ||2^. In our approach we use momentum dependent meson-quark 
interactions which provides for an effective cut off of the loop integration. We would also like 
to mention a recent investigation ||21[] where the same quark-meson Lagrangian employed in 
|T^ was used. The authors of [^| employed dipole vertex to describe various leptonic and 



semileptonic decays of both the heavy-light mesons and the i?c-nieson. 



In this paper we follow the strategy adopted in Refs. ||16|JT7|| . The basic assumption on 
the choice of the vertex function in the hadronic matrix elements is made after transition 
to the momentum space. We employ the impulse approximation in calculating these matrix 
elements which has been used widely in phenomenological DSE studies (see, e.g., Ref. [0). 
In the impulse approximation one assumes that the vertex functions depend only on the 
loop momentum flowing through the vertex. We present a general method which greatly 
facilitates the numerical evaluations that occur in the Feynman-type calculations involving 
quark loops (see also [0,^]). 

The basic emphasis of this work is to study leptonic and semileptonic decays of the B^. 
meson. We use Gaussian vertex functions with size parameters for heavy-light mesons as 



in Ref. IjT^- In this paper we limit our attention to the basic semileptonic decay modes of 
the i?c-nieson. A new feature of our calculation is that we also discuss semileptonic decays 
involving the r-lepton. We discuss in some detail how our quark loop calculations reproduce 
the heavy quark limit relations between form factors at zero recoil. Explicit expressions for 
the reduced set of form factors in this limit are given. 



II. MODEL 



We employ an approach |jT2[ based on the effective interaction Lagrangian which describes 



the coupling between hadrons and their constituent quarks. For example, the coupling of 
the meson H into its constituents qi and q2 is given by the Lagrangian 

Cintix) = gHH{x) jdxijdx2^H{,x,Xi,X2)q{xi)VH\Hq{,X2) . (2.1) 

Here, \h and Th are Gell-Mann and Dirac matrices, respectively, which entail the flavor 
and spin quantum numbers of the meson H. The function is related to the scalar part 
of Bethe-Salpeter amplitude and characterizes the finite size of the meson, is invariant 
under the translation $i/(x + a, xi + a,X2 + a) = $//(x,xi,X2) which is necessary for the 



Lorence invariance of the Lagrangian ( |2.1| ). For instance, the separable form 



$h(x, Xi, X2) =5{x- ^i±^) /((xi - X2f) (2.2) 

has been used in for pions with /(x^) being a Gaussian. The straightforward general- 
ization of the vertex function (2.2) to the case of an arbitrary pair of quarks with different 
masses is given by 
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^h{x,xi,X2) = 6[x ■ / xi-xsr . 2.3) 

\ mi+m2 J 



The authors of used a dipole form for the Fourier-transform of the function /(x^). Here 
we follow the slightly different strategy as proposed in Refs. []TB|,|T^ . The choice of the vertex 



function in the hadronic matrix elements is specified after transition to momentum space. 
We employ the impulse approximation in calculating the one-loop transition amplitudes. In 
the impulse approximation one assumes that the vertex functions depend only on the loop 
momentum flowing through the vertex. The impulse approximation has been used widely in 
phenomenological DSE studies (see, e.g., Ref. [|16|)- The final results of calculating a quark 
loop diagram depends on the choice of loop momentum flow. In the heavy quark transitions 
discussed in this paper the loop momentum flow is, however, fixed if one wants to reproduce 
the heavy quark symmetry results. 

To demonstrate our assumption, we consider the meson mass function defined by the 
diagram in Fig. |I|. We have 

Unix - y) = Jdxijdx2^ nix, xi,X2) jdyijdy2(^Hiy,yi,y2) (2.4) 

•tr {S{yi - Xi)ThS{x2 - y2)'^H} ■ 
Then we calculate the Fourier-transform of the meson mass function (12. 41"). 



flnip) = J e-'P^Unix) = (2.5) 

/(l^/(0^/(^'^^^~^' -^2)*^^^^' ^2)*^ {S{h)THSMTH} . 

The Fourier-transform of the function $(xi,...,x„) which invariant under the translation 

Xi ^ Xi + a can be written as 

71 

$(gi,...,g„) = dxi.. dxnC $(xi,...,Xn) (2.6) 

71 

■{2tt)^6 [^(lij ■ Jdxi.. JdxnS l^Xij e $(xi, ...,x„) 



^{27r)'5iY.qA^{q„...,q.^ 



n—l I 



\i=\ / 

Using this property one finds 

^h{v) = /^0|(fc,p)tr{5(^+ ^)THSif^)TH}. (2.7) 

Here, we assume that the vertex function (pH depends only on the loop momentum k. 
Besides, we assume that (pH is analytical function which decreases sufficiently fast in the 
Euclidean momentum space to render all loop diagrams UV finite. 
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The coupling constants gn is determined by the so called compositeness condition pro- 
posed in 1 13] and extensively used in [14|. The compositeness condition means that the 
renormalization constant of the meson field is equal to zero 



1 - 



0. 



(2.8) 



where 11^ is the derivative of the meson mass function defined by the diagram in Fig. |l| 



p2 



d^k 

An'^i 



(2.9) 



(2.10) 



For simplicity, we extract the factor l/47r^ from the definition of the meson mass operator. 
We use the local quark propagators 



1 



nii 



(2.11) 



where rrii is the constituent quark mass. As discussed in [|l^], we assume that trh < mq^+niq^ 
in order to avoid the appearance of imaginary parts in the physical amplitudes. This is a 
reliable approximation for the heavy pseudoscalar mesons, the above condition is not always 
met for heavy vector mesons. As discussed in Sec. VI we shall therefore employ equal masses 
for the heavy pseudoscalar and vector mesons in our matrix element calculations but use 
physical masses for the phase space. 



III. A METHOD FOR THE EVALUATION OF ONE-LOOP DIAGRAMS WITH 

ARBITRARY VERTEX FUNCTIONS 

For the present purposes one has to evaluate one-loop integrals of two- and three-point 
functions involving tensor integrands and product of vertex functions. In this section we 
describe a general method to efficiently enact these calculations for the general case of n- 
point one-loop functions. We note two simplifying features of our integration technique. 
The arising tensor integrals are reduced to simple invariant integrations. The sequence of 
integrations is arranged such that the product of vertex functions is kept to the very end 
and allowing for a full fiexibility in the choice of vertex functions. 

We consider a rank s tensor integral in the Minkowsky space as it appears in a general 
one fermion-loop calculation of a n-point function (see the diagram in Fig. One has 

^K"" = / (3.1) 

^ n [m1-{k + kY] 

i=l 

The outer momenta pj {j = l,...,n) are all taken to be incoming. The momenta of the 

i 

inner lines are given hj k + k with k = Pj such that = 0. The maximum degree of 
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the momentum tensor in the numerator arising from the n fermion propagators is Smax = n. 
We employ the impulse approximation dropping the dependence on the external momenta 
inside the vertex functions and denote the product of all vertex functions by J-'{—k^). 
Using the a-parameterization of Feynman one finds 

n 

where P = J2 caU-, and Dn{a) = OiiOj dij with dij = (1/2) [mj + m| — (/j — /j)^]. 

i=l ij 

Next we use the Cauchy integral representation for the function T{—k^) leading to 

The new denominator factor is then included again via Feynman parameterization giving 
= r(n + 1) / dpp--' J d-a6{l -ta.)J^f "^^^^'^^ 



27ri 



X 



- n+l ■ 

1 -/?) C - (A; + /?P)2 + -/?(!-/?) P2 

One then factors out the (1 — /3) in the denominator and shift the integration variable k to 
k' ^{k + (3P)/^/T^ to obtain 

(yr^k- ppf"... [-^T^k- i3py' 

^ [C - A;2 + zY^^ ■ 

The contour integral can be done again by Cauchy's theorem. On substitution oi (5 — 
tj (1 + one then has 

S'"'^^ - i-r J dtjY-^, J d-c.s{i - E «o / -?^^"^(-^' + ^) ' 

where ^("^ denotes the n-th derivative of the function T and where 

Ri'^^^k^ —Pi", z^tDJa) —P^. 

VT+t l + t ' ^ l + t 

The momentum integration of the tensor integral can be trivially done by invariant integra- 
tion. Finally we go to the Euclidean space by rotating kg ik^ which gives k^ — > —k^ = u, 
then one encounters the scalar integrals 
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j.n— 1 



:i+t) 

The u-integration can be performed by partial integration and one finally obtains 



oo 1 



/J.IL J. /I /(. 

(TT^ y '^""'^^^ ~ ^ «.)^(— 2)(z) . (3.4) 

Since mmax = [^^72] Eq. ( p.4|) holds true for all n and m except for the case n = 2 and m = 1. 
In this case we have 

/[2,i] = 2 y" c^^ ^ X3 y c?^a5(l - / duJ^{u) 



1=1 



oo 2 2 

= / (iTt) / '^'"'^^^ " ? «.)zj^(^) . (3.5) 



where = dz{t)/dt. 

One has to remark that the integration over the a— parameters in Eq. ( |3.4| ) can be done 
analytically up to a remaining one-fold integral. However, the ease with which the numerical 
a-integrations can be done does not warrant the effort of further analytical integrations. 
Using the integration techniques described in this section all necessary numerical integrations 
encountered in this investigation can be performed within minutes using a fast modern PC's. 



IV. HADRONIC MATRIX ELEMENTS 



A. Quark-meson coupling constants 



As already discussed in Sec. II, the quark-meson coupling constants are determined by 
the compositeness condition Eq. ( p. 81) . The derivatives of the meson-mass functions can be 
written as 



/^np(/ 



1 



d rd'^k 



2p2^ dp'^JAnH^^ 
d^k 



7^53(^)7'5i(^+ ^) 



2p^ J AttH^^^ ' 



-i''S^{}i)-i''S^{}t+ /) ^S^ 



(4.1) 



_d_ 

dp'^ 



1 1 



d 



p 



3 dp 



p^p^ 

p2 



d^k 



0y(-A:')tr 



(4.2) 



p2 



d'^k 2 

2p2 J 4^*^^ 



-A;2)tr 
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The evaluation of the integrals is done by using the method outlined in Sec. III. The 
compositeness condition reads 



^iV^ = l, where Nh = ^IIh{v%- 



(4.3) 



H 



1 °° t 2 ^ 

- dti I da a 

2J \l + t) J 





H 



-V 



1 



1 + t 

at ^2 



4-3- 



1 + t 



2 - 



at 
at 



T'p{z) \ 2mim3 + 



at 



2 2 

— 2mim3 + p 



1 + t 



-V 



1 



at 



1+t 
at ^2 



1 + t 



1 + t 
at 



T'y{z) \ 2mim3 + 



at 
TTt 



— 217111713 + 



1 + t 



Here, mi stands for the heavy quark {b or c) and m2 for the light quarks {u, d, s) in the case 
of heavy-light systems, and for c in the case of double-heavy systems. The function J-'h{z) 
is the product of two vertex functions J^h{z) = (P'h{z) with 



z = t (ami + (1 — a)m'l — a{l — a)p^^ — — p^. 

1 + t 

B. Leptonic and radiative decays 

The matrix elements of the leptonic and radiative decays are defined by the diagrams in 
Figs. H-|3 and given by 



-Cv 



rd^k 






M- 




"d^k 


Any 





fpp' 



(4.4) 



0y(-P)tr 



-mvCvfve*' 



(4.5) 
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WP-n ^ €2 Qi q2 , 



(4.6) 



^9v9p 


rd^k 


47r2 J 





l'S,{}t+i^)fSMf,S,{}i+i) 



(4.7) 



For ease of presentation, the expression for the V P7-decay is given for neutral-flavored 
mesons. Using the integration techniques described in Sec. Ill one then arrives at the follow- 
ing analytical representation of the various one-loop matrix elements 



^ '0 



+ (mi - ma) 



at 



1 + t 



(4.8) 



ami + (1 ~ a)m^ — a{l — a)p 1 — - — 

1 "h ^ 



p 



mv 47r2 J (1 + ty 
^ ^0 



2 1 , at 



1 - 



at 

iT7 



(4.9) 



zy = t (rrig - a{l - Q;)p^) - j 
_L . 47r 



r(V-^e+e-) = ^ — /^C^ Cy = el {V = <!> , J/^ ,T). 
my ^ 



gpjj = Cp^^ ' ' ^ / ^^(iTt) / ^^"^(1 ~ ^i) (-</'p(^o)) , (4.10) 



2:0 = t ^m^ - Q;ia2P^) - ^ ^ ^ ^1^2^^ , 



t (r 



TT 



r(p ^ 77) = -a TJipgp^^ , 



(7 = 



gvp-i = C'yp^ • mq • y ^^(7+7) / ^^"^(^ ~ ^'^^) (~-^Vp(^v-p)) , (4.11) 
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zvp = t (niq — aitts my — aia2 rrip^ — («i(ai + cx2)my — aia2^p)) j 

Pl) = Y^ml (l - ^) \Ip, Cj,^,^, = 2e, . 



The electric quark charges Cq are given in units of e. 

C. Semileptonic form factors 

The semileptonic decays of the i?c-nieson can be induced by either a beauty quark or 
a charm quark transition. In the relativistic quark model, the hadronic matrix element 
corresponding to b-decay is defined by the diagram in Fig. ^ and is given by 



^9p9h 


rd^k 


47r2 J 





MtiPip) - H{p')) = h-n:FpH{-e)ir 



, (4.12) 



where J^ph = (pp ■ 4>h, ^p = 1^, and Ty = —i with e* ■ p' = 0. For the b-decay case one 
has the CKM-enhanced decays 

b —>■ c : {rjc, J / ip) I^p mi = nib , m2 = m^ = mc , 

and the CKM-suppressed decays 

b ^ u : (D'^, D*'^) l^u mi = mb , m2 = mu , m^ = mc- 

The c-decay option of the i?c-meson is represented by the Feynman diagram in Fig. |^ 
which gives 



Mi!{P{p) ^ H{p')) = J^J—^^J^pH(-k')tr^^'SsmO'S2i^+ ^)THSii^+ . (4.13) 

Again one has the CKM-enhanced decays 

c — >■ s : B'^ — i> {B^, B*^) l^v mi = mb , m2 = m^ , m^ = mc- 
and the CKM-suppressed decays 

c ^ d : B^ {B^i B*^) l^v mi = mb , m2 = ma , m^ = mc ■ 

It is convenient to present the results of the matrix element evaluations in terms of invariant 
form factors. A standard decomposition of the transition matrix elements into invariant 
form factors is given by 

M^{P{p)^P'{p')) = U{q') {p + p'r + f-{q') (p-pT (4.14) 

and 
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iMf'iPip) V{p')) = -g^"'e*''{mp + mv)Ai{q^) + {p + p'Y p ■ e 



nip + iTLv 



(4.15) 



+ {p-p'Yp-e* 



mp + my mp + mv 



The various invariant form factors can be extracted from the one-loop expressions (|4.12|) 
and ([4.131) by using the techniques described in Sec. III. One finds that the form factor 
integrands factorize into a common piece times a piece specific to the different form factors. 
One can thus write 



a, 



(4.16) 



where F = f±,Ai,V. For the ^0 b c,u form factors /+ one has 

t 



^ 1 
= Tpp{zh) 

f+ 



l + t 



4 - 3(ai + 02) 



l + t 



- J^pp{zb)^{mi + m2)m3 



Y^^(^-(q;i + a2){mim3 + 7712-023 - 7771777,2) + aip"^ + 0^2 p'^^ 

2 - (ai + 0:2) Y~j~l) (^"^ ~^ "2) («! + «2 p'^) - aia2 | 



l + t 



\i=i 



aim^ — aia^ip^ — 020:3^'^ — 010:20'^ 



l + t 



P"^ , Pb = aip + a2p' ■ 



For the corresponding c — >■ s, d form factor one has 



- -^Pp{Zc)-—— 

/+ l + t 



1 + 3ai 



l + t. 



+ J'pp(^c)|"^2m3 

+ Y^t ('^l("^l"^2 + 777x7773 - 7772^^3) + «2 + (j^^^ {^iV^ " «2 9^)) 

- (^Y^^) '^1(^0^1(0^1 + "2)p^ - aia2p''^ + Q;2(ai + ^2) , 

= t ai777,^ - aia3 - 02^3 - 0L\012 j - -^c^ ' = OL\V + "2? • 

Expressions for the remaining 0~ 0~ and 0~ ^ 1~ form factors /_, Ai and V are given 
in the Appendix. The masses 777, (i = 1,2,3) appearing in the form factor expressions are 
constituent quark masses with a labelhng according to Eqs. ([4.12|) and ( |4.13| ). The values 
of the constituent quark masses as well as the vertex functions entering the form factor 
expressions will be specified in Sec. VI. 
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For the calculation of physical quantities it is more convenient to use helicity amplitudes. 
They are linearly related to the invariant form factors [2^. For the 0~ — 0~ transitions one 
has 



i^o(a = ^A(g^), (4.17) 

H,{q') = -^{(m^ - ml,) U{q^) + /_(g2)} . (4.18) 



For the 0^ 1^ transitions one has 



H^q') = -{mp + my)A,{q') T 7 V{q') , (4.19) 

[mp + mv) 

Ho{q') = ^ -{ml -ml- q'){mp + my)A,{q') + f A,{q') , (4.20) 

2mv\/q^y mp + my ) 

HM^) = ^^^|-(mp + mv) A,{q') + {mp - my) ^^(g^) + A,{q') ] , (4.21) 

mvVQ t mp + my J 

^X{ml,ml,q^) _ [(^2 _ g2)(^2 _ ^2)]i/2 



^vVg mp + mv 

where 

P 

2mp 2m p 

with = {mp ± mu)"^. 

Then the partial helicity rates are defined as 

d?-(2^'^^^'' ■ 12m|,g2 ' )l ' ^-±,0,t, (4.22) 

where Vff is the relevant element of the CKM matrix, mi is the mass of charged lepton. 
Finally, the total partial rates including lepton mass effects can be written as 



dq^ ~^^ + 2g2^ dq^ ^ ^ 2q^ dq^ ' ^ 



dr 



PV 



'1 + ^) 
^ 2q^' 



dr 



PV 



dg2 



+ 



dT 



PV 



+ 



drr 



dq^ 



ml dvr 



+ 3-^ 



2g2 dg2 



(4.24) 



In the following we shall present numerical results of the total decay widths, polarization 
ratio and forward-backward asymmetry. The relevant expressions are given by 



, 2 



a 



To 



r, +r_ 



A 



FB — -7 



3 r_ 
4 



(4.25) 
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V. HEAVY QUARK SPIN SYMMETRY 



Our model allows us to evaluate form factors directly from Eq. ( [4.16D without any ap- 
proximation. However, it would be interesting to explore whether the heavy quark spin 
symmetry relations derived in Ref. can be reproduced in our approach. As was shown 
(see, for instance, IIT^) our model exhibits all consequences of the spin-flavor symmetry for 



the heavy-light systems in the heavy quark limit. For example, the quark-meson coupling 
and leptonic decay constants behave as 

2 oo , 



3N, 



H 



fn-^^- J-^ [du[V^ - E]Mu - 2EV^) ri^'^l^^ , (5.2) 
^/ml \ 27T^Nh ml + M — 2_£/a/m 

in the heavy quark limit: = m'jj = (mi + when mi 00. Eqs. ( p.l| ) and ( ^.2|) 
make the heavy quark mass dependence of the coupling factors gn and fn explicit since we 
have factorized the coupling factor contributions into a heavy mass dependent piece and a 
remaining heavy mass independent piece. Moreover, Eqs. (|5.1|) and ( |5.2|) show gn and fn 
scale as ■m\^'^ and m^^^'^, respectively. 

As is well known (see Ref. |Q), heavy flavor symmetry cannot be used for hadrons 
containing two heavy quarks. But one can still derive relations near zero recoil by using 
heavy quark spin symmetry. 

First, we consider the semileptonic decays Be — > Bs{B^)e^v and B^ B*{B*^)e'^u 
which correspond to c-decay into light s and d quark, respectively. Since the energy released 
in such decays is much less than the mass of the b-quark the four- velocity of the 5c-nieson 
is almost unaffected. Then the initial and final meson momenta can be written as 

p = rriB^v p' = mBV + r 

where r is a small residual momentum [v-r = — r^/ {21713))- The heavy quark spin symmetry 
can be realized in the following way. We split the B-meson masses into the sum of b-quark 
mass and binding energy 

rriB^ = mp = mi + Ei, tub = run = rrii + E2. 

Then we go to the heavy quark mass limit mf, = mi ^ 00 in which the b-quark propagator 
acquires the form 

' ^ '^'^ (5.3) 



mi- ^- ^ ~2{kv + El)' 

The decoupling of the c-quark spin allows us to reliably neglect the fc-integration because k 
is small compare to the heavy c-quark mass. One has 

mg— f m3 
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As a consequence, the hadronic matrix element describing the weak c-quark decay simphfies: 
J2m£^2ml 1 rd'k , ,^,, tr[0^(m2+ ^+ i)VHi\^ 



where q = p — p' = {mp — mH)v — r={{Ei — E2)v — r = AEv — r and m2 stands for the 
hght quark mass (m2 = mg or rrid). One has to emphasize that all above approximations 
are valid only close to the zero-recoil point q"^^^ = AE"^. Recalling the transversality of the 
final vector meson field p' ■ e*=mBV ■ +r ■ e* = and applying the integrations as described 
in Sec. Ill, one finds 



r T f at-AE\ „ 1 „ 



t{...}pv = -le^^'-^e* ^"r^YTt V"' ^ j - ITt '''^ ' ^ ' 

Here, Zc = {at^ /{l+t)){a + 2AE) +t{ml-2aEi) - {t/{l + t)) AE"^. It is readily seen that 
the amplitudes of c-decay in the heavy quark limit are expressed through two independent 
functions 

r T f at-AE\ , . 1 

To complete the description of the heavy quark limit in the c-decay modes, we give the 
expressions for the form factors in this limit. One has 



F(,U ^^^^ . I (5.7) 



Np-Nh 







where F = f±,Ai,V. The form factor specific pieces are given by 

r T 1 / at \ , , 11 



2mim3 V 1 + tJ 1 + t 

1 If at-AE 

= -- — — — ■ — ■ [^2 + 



mp + my \ 1 + t 

. , . , mp + mv 1 1 



{...} 



mp + my 



As 



4(^-1) 



1 + t ms 



Superficially it appears that the form factors /+ and Ai are suppressed by a factor of 1/mi. 
However, they must be kept in the full amplitude to obtain the correct result in Eq. ( p. 61) , 
for instance, one has 
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UiP + pT + f4p- P'Y = (2mi/+ + ^Ef^)v^ + (/+ - /_)r^ 

A similar analysis applies to the h u decays B,. —>■ {D^, D*^) e^u. Again the heavy 
quark symmetry analysis is only reliable close to zero recoil where the u-quark from the 
h ^ u decay has small momentum. One has 

V2m732m^ 1 [d^k , ,^,, tr[7^rH(m2+ ^+ ^')0^^{l+ ^)] 



where q = p — p' = (mi + Ei — mH)v — r. The light quark mass m2 in Eq. ( ^.8|) is the 
u-quark mass. One finds 



oo , oo 



MS = . ± . / « / <i„^,„(,,){...},„ (5.9) 



with 



rriH — <yt\ 1 



l + t 



i{..\py = -ie^'-»(^e*•'v"r'^ ( ms + — — —) e""''^ — v^" e* ■ r . 

V l + t J l + t 

Here, Zb = [at^ / (1 + 1)) (a + 2mj^) + 1 (rn\ — 2aEi) — {t/ (1 + 1)) mj^. Again, the amplitudes 
for the h ^ u decays are expressed through two independent functions. The expressions for 
the form factors in the heavy quark limit close to zero recoil read 

2 ^ v^2mp ■ 2mH °f dtt 



^ / /I , / daJ^PH{zc){...}F , (5.10) 
Np ■ Nh ^ ' 



where E = f±, Ai,V and where 

1 1 



{•••}/+ = -{•••}/- 



2m3 l+t 

rriH — at 



= ■ [7712 + 

771 p + 77lv 771^ 

{■■■]a2 = -{■■■]a3 = {■■■}v 



77lp + 77lv 771^ V l+t 
771 p + 77lv 



rriim^ l + t 



Note that one needs to keep the next-to- leading term in the sum (/+ + /-) to obtain the 
above amplitudes. 

The hadronic matrix elements of the 6 — > c decays B^. —>■ {7]c, J/tp)e'^i' simplify signifi- 
cantly in the heavy quark limit. In this case both b- and c-propagators may be replaced by 
their heavy quark limit forms in Eq. ( [5.3| ) with the same velocity v. Again, the results will 
be valid only near zero recoil. One has 
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ms y AttH ' ' [-2A;t; - 2Ei] [-2kv - 2E2 



'Np ■ Nh 

where p = (mi + Ei)v, p' = {m^ + E2)v + r. One finds 

00 1 



^^^^^ y^^^ 1 _ r I ^^jr^^f_2^^aE, + {l-a)E2)){...}pH (5.12) 



{...}PP = +2^;^ ^{...}py = -2e*r 

The form factors are written down 

00 1 

V2mp ■ zniH 



F(gLx) ^ ^^-^2^=^ ■ [du [ daJ^PH(u-2V^{aEi + {l~a)E2)){...}F (5.13) 





where F = f±,Ai,V. We have 

_ mi + m3 r I _ ^1 -^3 

^•••^^+ ~ 4mimi ' ~ Am^ml ' 

1 1 r I _ r I _r I _ mp + my 

Thus, our quark loop calculations reproduce the heavy quark limit relations between form 
factors obtained in |^ near zero recoil. Moreover, we give explicit expressions for the reduced 
set of form factors in this limit. 



VI. RESULTS AND DISCUSSION 

Before presenting our numerical results we need to specify our values for the constituent 
quark masses and shapes of the vertex functions. As concerns the vertex functions, we found 
a good description of various physical quantities |T^ adopting a Gaussian form for them. 
Here we apply the same procedure using (pnik"^) = exp{— /c^/A|^} in the Euclidean region. 
The magnitude of characterizes the size of the vertex function and is an adjustable 
parameter in our model. We reiterate that all the analytical results presented in Sec.V are 
valid for any choice of form factor (pnik^)- For example, we have reproduced the results of 



21| where dipole form factor was adopted by using our general formula. 

In [0 we have studied various decay modes of the vr, K, D, Dg, B and Bg mesons. 
The A-parameters and the constituent quark masses were determined by a least-squares fit 
to experimental data and lattice determinations. The obtained values for the charm and 
bottom quarks (see Eq. ( |6.1| ) allow us to consider the low- lying charmonium {rjc and J/ip) 
and bottonium (T) states, and also the new-observed i?c-nieson. 

(6.1) 



0.235 0.333 1.67 5.06 
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Basically we use either the available experimental values or the values of lattice simula- 
tions for the leptonic decay constants to adjust the size parameters Ah- The value of fs^ 
is unknown and theoretical predictions for it lie within the 300-600 MeV range. We choose 
the value of /^^ = 360 MeV, being the average QCD sum rule predictions, for fitting A^^. 
The obtained values of Ah are listed in Eq. (p.2|) as well as the values of fn in Table I. 



Arr Ak Ap Ap^ Aj/^ Ab Ab, Ab, At , . 

1.16 1.82 1.87 1.95 2.12 2.16 2.27 2.43 4.425 ^ ' 

The values of Ah are such that A^. < A^ if < Tnj. This corresponds to the ordering 
law for sizes of bound heavy-light states. 

The situation with the determination of A^^ is quite unusual. Naively one expects that 
A^^ should be the same as Ajj^. However, in this case the value of the r]c —>■ 77 decay width 
comes out to be 2.5 less than the experimental average. The experimental average can be 
reached only for relatively large value of A^^ = 4.51 GeV. Note that the values of the other 
observables {J/ip — »■ rjc'J and — »• rjclu decay rates) are not so sensitive to the choice of 
A„ : 

Br(r/c ^77)= 0.031 (0.012) % , expt. = (0.031 ± 0.012) % , 
Br( J/V^ ^ r/c7) = 0.90 (1.00) % , expt. = (1.3 ± 0.4) % , 

Br(5e ^ r]Ju) = 0.98 (1.02) % . 

The values in parenthesis correspond to the case of equal sizes for the charmonium states. 

We concentrate our study on the semileptonic decays of the i?c-nieson. To extend the 
number of modes, we consider also the decays into the vector mesons D*, B* and B*. We 
will use the masses and sizes of their pseudoscalar partners for the numerical evaluation of 
the form factors to avoid the appearance of imaginary parts in the amplitudes. Such an 
assumption is justified by the small differences of their physical masses. 

In Figs, (pfjl^) we show the calculated dependence in the full physical regions of 
the semileptonic form factors of the CKM-enhanced transitions B^. rjc, B^ —>■ J/ip and 
Be Bg, Be ^ B*. The values of form factors at maximum and zero recoil are listed in 
Tables II-IV. The comparison of the exact values of form factors at zero recoil and those 
obtained in the heavy quark limit is given in Table V. Our results indicate that the corrections 
to the heavy quark limit at the zero recoil point = q^^^ can be as large as a factor of two 
in 6 — c transitions and a factor of almost five in b — u and c — d transitions. This is not 
so surprising considering the semileptonic decays of the D meson where similar corrections 
can amount to a factor of two |[TB[ . 

The form factors can be approximated by the form 

/(^') = 1 2, 2^^°i ,2/ 2,2 (6-3) 

1 - g2/m4 - 6 ■ {qymi^y 

with the dimensionless values of /(O) given in Tables II and III. Note that the form factor 
f+{q^) for the B^. —>■ rj^ transition rises with as is appropriate in the time-like region. 
When plotted against uj = p ■ p'/ {rrii^ ■ mout) = ("^fn + ^lut ~ 'f) I (2 "^m ■ "^out), they would 
fall with UJ as one is familiar with heavy quark effective theory. 
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It is interesting that the obtained values of m^^ for the CKM-enhanced transitions (see 
Table VIII) are very close to the values of the appropriate lower-lying {qq') vector mesons 
{niB* ~ rriB^ = 6.4 GeV for b ^ c, mj:,* = 2.11 GeV for c s.). The parameter S charac- 
terizes the admixture of a g^-term in the denominator. Its magnitude is relatively small for 
all form factors of the CKM-enhanced transitions except Ai for Be ^ J/ip transition and 
A2 for Be —>■ B* which have a rather flat behavior. This means that those form factors can 
be reliably approximated by a vector dominance form. However, one cannot approximate 
the form factors for the CKM-suppressed transitions by a pole-like function only. 

We use the calculated form factors in Eq. (|4.25|) to evaluate the branching ratios for 



various semileptonic B^ decay modes including their r-modes when they are kinematically 
accessible. We report the calculated values of a wide range of branching ratios in Table VI. 
The results of other approaches are also given for comparison. The values of branching ratios 
of the CKM-enhanced modes with an electron in the final state are of order 1-2 %. The 
values of branching ratios of the CKM-suppressed modes are considerably less. The modes 
with a T-lepton in the final state are suppressed due to the reduced phase space in these 
modes. To complete our predictions for the physical observables we give in Table VII the 
values of the polarization ratio and forward-backward asymmetry for the prominent decay 
modes. 
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VII. APPENDIX 



In this Appendix we list the remaining form factor expressions appearing in the curly 
brackets in Eq. (|4.16| ) which have not been listed in the main text, 
b-decay: 



-J^pp{Zb) 



1 + t 



3(ai - ^2) 



1 + t 



+ J^'pp{zb)\ {rrii - 7x12) mg 



l + t 



(tti - a2)("^i"^3 + m2m3 - mim2) + aip^ - a2p' 



12 



-(ai - 02 



[ai + a2){aip + a2P )-aia2q 



l + t 



mp + my L 



PV 



l + t 



(mi + 2m2 - ma) - J^'pyizb) 



mim2m3 + ^(p^ - g^)m3 
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1 t 
+- 



; f (aimi + {2ai + 0:2) m2 - (Scki + 02)1713) 



21 + t 

+ ((cKi + 2a2)mi + 021712 — (cKi + 3a2) ma) p'^ 
- (cKimi + Q;2m2 - (cti + q;2) ma) g^^ 

-(mi + m2 - m3)^(Q;i + q;2)(q;iP^ + 0:2 p'^) - Q;iq;2 (^j 



+ t 



A2 



{mp + mv) 



^Pvi^b) 



-ma - 



1 + t 

2 



aimi + Q;2m2 — (3q;i + q;2) ma 



+2(mi - ma) cti (ai + 0:2) (y^) } 



A3 



{mp + mv) 



^pvi^b) 



ms + 



t 



+2(mi - ma) ai {ai - 0:2) 



l + t 

t -2 



aimi + Q;2m2 + {ai — a2) ma 



l + t 



(mp + my) 



-^pvi^b) 



ma + 



ai(mi - ma) + 02 (m2 - ma) 



We use the abbreviations 



Pb = otip + a2p' , Zb^t ~ CKiCCap^ - Q;2Q;ap'^ - aia2q^ - Ph 



c-decay: 



" 1 

= 2>Tpp{zc)-— 
f- l + t 



t 



1 - (ai + 2q;2) 



l + t. 



+ T'pp{zc)^-2mim3 + m2ma 



l + t 



(^{ai + 2a2){mim2 + mima - m2ma) - 2{ai + q;2)p^ + 2a2p'^ - a2q 
+ (^Y^T^) (("^"^i + 6q;iq;2 + 2q;2) p^ — 2a2{oL\ + q;2) p'^ + q;2(2q;i + 3q;2) (f^ 
(«! + 2ci;2) {oLx{ax + 0:2) - Q!iQ;2P'^ + ol2{olx + 0:2) g^^. 
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2 f 1 



mp + my I 
1 t 



1 + r 



-mim2m3 + + - p'^) ma 



(ttimi - {2ai + a2) m2 - (3q;i + a2)^3) P 



^21 + i 

+ (-aimi + Q;2m2 + (ai + 0:2) ms) p'^ 
+ ((cKi + 2Q;2)mi - Q;2m2 - (ai + 3q;2) ms) 

- ^Y^-^^ ("^1 -^2- ms) ^Q;i(ai + 0:2)]?^ + 02(0:1 + 02) - q;iQ;2p'^^|- 



|...|^ =(mp + my) -T'py{zc) jma + 



t 



2 



ttimi — a2m2 — (3q;i + a2) ma 



-2(mi - ma) ai (ai + q;2) (j^) } 



{mp + my) 



^'pvi^c) 



-Sma + 



1 + t 



— (3q;i + Aa2)mi — a2m2 + (Scti + 7q;2) ma 



+2(mi - ma) (ai + 2^2) ("i + ^2) 



/ t 



(mp + my) 



"•^py(^c) 



ma + 



1 + t 



Q;i(mi - ma) + Q;2(m2 - ma) 



Here we have used the abbreviations 



Pc = Oiip + a2q , Zc^t aim^ - aia^p^ - a2as - aia2p''^ - Pc- 
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TABLES 

TABLE I. Leptonic decay constants fn (MeV) used in the least-squares fit. 



Meson 



This model 



Other 



Ref. 



TT 



+ 



131 



130.7 ±0.1 ±0.36 



Expt. 



160 



159.8 ± 1.4 ±0.44 
^ 



Expt. 



191 



191] 



28 



192 ± 11+8*LV^ 
l94:t\o ± 10 



-28 



Lattice |2|j24| 
Lattice |25tl 
Lattice |26 



206 



206 



210 ±9 



+25+17 
^9-1 



213l|i ± 11 



=7r 



Lattice 
Lattice 
Lattice 



2|,|24| 



2|,|24| 
2| 



172 



172; 



-31 

157 ± ir 

164+lf ± 8 



'25+23 
9-0 



171 ± lOttT 
185^1^ ± 9 



Lattice 
Lattice 
Lattice 



196 



Lattice 
Lattice 



Br 



479 
500 
512 
687 
480 
432 

400± 20 
300 

360± 60 
300± 65 
385± 25 
420(13) 



logarithmic potential [27| 
Buchmiiller-Tye potential |27] 
power-law potential 
Cornell potential |^ 
Potential model 
QCD-inspired QM [ 
QCD spectral SR l2 
QCD SR 13 



QCD SR H 



QCD SR |32 
QCD SR 13 



Lattice NRQCD m 



Br 



360 



360 



Our average of QCD SR 



404 



405± 17 



Expt. 



T 



711 



710±37 



Expt. 
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TABLE II. Predictions for the form factors at = and = q'^gx foi" ^ P decays. 



p 


/+(0) 


/ + (^max) 


/-(O) 


/— (^max) 




0.76 


1.07 


-0.38 


-0.55 


D 


0.69 


2.20 


-0.64 


-2.14 


B 


-0.58 


-0.96 


2.14 


2.98 


Bs 


-0.61 


-0.92 


1.83 


2.35 



TABLE III. Predictions for the form factors at = for Be ^ V decays. 



V 


Ai(0) 


^2(0) 


^3(0) 


F(0) 


J/^P 


0.68 


0.66 


-1.13 


0.96 


D* 


0.56 


0.64 


-1.17 


0.98 


B* 


-0.27 


0.60 


10.8 


3.27 


B*s 


-0.33 


0.40 


10.4 


3.25 



TABLE IV. Predictions for the form factors at = g^ax B^—^V decays. 



V 


^l(5max) 




^3(gmax) 






0.86 


0.97 


-1.71 


1.45 


ir 


0.85 


1.7(3 


-3.69 


3.2(3 


B* 


-0.42 


0.49 


18.0 


5.32 




-0.49 


0.21 


15.9 


1.91 



TABLE V. Comparison of the form factors at the zero recoil point = g^gx calculated in the 
heavy quark limit with exact results. 



H 


/+ 


/- 


Ai 


A2 


^3 


V 




1.07 


-0.55 


0.86 


0.97 


-1.71 


1.45 


ijc,J/i^ (HQL) 


0.70 


-0.35 


0.37 


0.69 


-0.69 


0.69 


D, D* 


2.20 


-2.14 


0.85 


1.76 


-3.69 


3.26 


D, D* (HQL) 


0.59 


-0.59 


0.18 


1.50 


-1.50 


1.50 


B, B* 


-0.96 


2.98 


-0.42 


0.49 


18.0 


5.32 


B, B* (HQL) 


-0.47 


1.67 


-0.25 


1.91 


21.47 


1.91 



25 



TABLE VI. Branching ratios BR(%) for the semileptonic decays Hl~^u, calculated with 

the CDF central value r{Bc) = 0.46 ps H. 



H 


This model 


1,0 





i 





i 


121] 




0.98 


O.SibO.l 


0.78 


1.0 


0.15(0.5) 


0.6 


0.52 




0.27 














J/ip eh' 


2.30 


2.1±0.4 


2.11 


2.4 


1.5(3.3) 


1.2 


1.47 


J ftp T V 


0.59 














L) e v 


U.Uio 




U.UUo 


U.UUb 


U.UUUo(U.UUz j 








0.0094 
















0.034 
0.019 




0.013 


0.019 


0.008(0.03) 






B^eu 


0.15 




0.08 


0.16 


0.06(0.07) 






B*^ e V 


0.16 




0.25 


0.23 


0.19(0.22) 






B^ev 


2.00 


4.0 


1.0 


1.86 


0.8(0.9) 


1.0 


0.94 


Bf e V 


2.6 


5.0 


3.52 


3.07 


2.3(2.5) 




1.44 



TABLE VII. The polarization ratio a and forward-backward asymmetry ApB- 



H 


a 


Afb 




1.15 


-0.21 




0.10 


-0.46 


B*° 


0.94 


0.35 




1.09 


0.29 



TABLE VIII. The numerical values of rrig^. (GeV^) and 6 in the form factor parameterization 
Eq. (P). 





U 


f- 


Ai 


A2 


^3 


V 


Be r]c, J/ip 




(6.37)2 


(6.22)2 


(8.20)2 


(5.91)2 


(5.67)2 


(5.65)2 




5 


0.087 


0.060 


1.40 


0.052 


-0.004 


0.0013 


Be Bg, B* 


^It 


(1.73)2 


(2.21)2 


(1.86)2 


(3.44)2 


(1.73)2 


(1.76)2 




5 


-0.09 


0.07 


0.13 


-107 


-0.09 


-0.052 



26 



FIGURES 



k+p 




FIG. 1. One-loop self-energy type diagram needed for the evaluation of the compositeness 
condition. 




Pn ^ \ P3 



I, = tp,, l„ = . 

J = l " 

FIG. 2. One-loop diagram with n-legs and arbitrary vertex functions. 



k+P 




FIG. 3. Quark model diagram for leptonic meson decays. 
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k 

FIG. 7. Quark model diagram for the semileptonic i?c-decays involving b ^ c,u transitions. 
Lower leg in loop is c-quark. 



k+P 




FIG. 8. Quark model diagram for the semileptonic i?c-decays involving c ^ s,d transitions. 
Upper leg in loop is b-quark. 
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1.1 





FIG. 9. g'^-dependence of the Be 
f-{q^) form factor. 



q^(GeV^) 

r7c form factors. Note that we plot the negative of the 
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